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ABSTRACT 

Two thirds of the F star members of the 12 Myr old j3 Pictoris Moving Group (BPMG) show 
significant excess emission in the mid-infrared, several million years after the expected dis- 
persal of the protoplanetary disc. Theoretical models of planet formation suggest that this 
peak in the mid-infrared emission could be due to the formation of Pluto-sized bodies in the 
disc, which ignite the collisional cascade and enhance the production of small dust. Here we 
present resolved mid-infrared imaging of the disc of HD 19 1089 (F5V in the BPMG) and con- 
sider its implications for the state of planet formation in this system. HD191089 was observed 
at 18.3|xm using T-ReCS on Gemini South and the images were compared to models of the 
disc to constrain the radial distribution of the dust. The emission observed at 18.3[J.m is shown 
to be significantly extended beyond the PSF at a position angle of 80°. This is the first time 
dust emission has been resolved around HD191089. Modelling indicates that the emission 
arises from a dust belt from 28-90 AU, inclined at 35° from edge on with very little emission 
from the inner 28AU of the disc, indicating the presence of an inner cavity. The steep slope 
of the inner edge is more consistent with truncation by a planet than with ongoing stirring. A 
tentative brightness asymmetry -pj = 0.80 ± 0.12 (1.8er) between the two sides of the disc 
could be evidence for perturbations from a massive body on an eccentric orbit in the system. 



1 INTRODUCTION 



The majorit y (^80 per cent) o f stars are born with protoplan- 
etary discs dHaischet all 1200 ll) which are massive and contain 
large amounts of gas, with an estimated gas to dust ratio of 100:1 
i Andrews & Williams 2005). These discs have a lifetime of -~6Myr 
l Haisch et al.ll200lh and at ages of a few Myr undergo a tra nsition, 



the nature of which is still unclear dMuzerolle et al. 2010l) . Some 
protoplanetary discs have inner holes dLuhman et aL 20 lOh but the 
systems retain an optically thick outer disc and are often called tran- 
sition discs to signify that this is a phase that discs go through. At 
later times (~10Myr) the remaining optically thick material in the 
disc has also disappeared, leaving optically thin debris discs. 

Debri s discs surround ~15 per cent of main sequence stars 
fsee lWvattl 12008) for a recent review). The dust particles in debris 
discs are generally small, and so are subject to radiation forces such 
as radiation pressure and Poynting-Rob ertson (P-R) drag from the 
central star ( Ba ckman & Par esce 1993) which remove grains from 
the system on timescales much shorter than the stellar lifetime so 
the dust must be continuously replenished. The collisional destruc- 
tion or sublimation of larger bodies (planetesimals) can create such 
dust particles. Effectively, the small dust can be used as a tracer 
for the presence of bodies many hundreds of kilometers in size that 
are otherwise undetectable. Inner cavities have also been observed 
in debris discs, either directly through images or indirectly through 
an absence of hot dust. The holes imply a lack of planetesimals in 
the inner region, as rapid collisional processing means that n one of 
the d ust created in the outer disc reaches the inner region dWvattl 
l2006h . 

Since gas giants must form in the gas-rich protoplanetary disc 



then it is possible that the planetary system is fully formed by the 
end of the transition phase. Imaging of debris discs at ages of ~ 
lOMyr shows evidence that planet formation in these systems has 
already finished and formed Jupiter mass planets. The ex istence 
of such planets was infe rred from debris disc observations ( Ouillen 
l2006tlWvatt et al . 2005) and subsequently confirmed through direct 
imaging of planets in systems such as the ~8M, U P planet imaged 



aroun d (3 Pic, a lOMyr old A star with a debris disc ( Lagran ge et al 
2009) and a <3M,„ P plan et imaged around Fomalhaut ( Kal as et al 



2008; C hiang et al.ll2009l) . 

However models indicate that the growth of solid objects in a 
debri s disc can continue for several 100 Myr dKenvon & Bromlevl 
1 20041) . In such self-stirred models a bright ring appears in an ex- 
tended planetisimal disc where Pluto sized bodies have just formed. 
These bodies stir the planetesimal belt into collisional destruction 
and increase dust production, causing the disc to brighten at this 
location. There is some evidence in favour of this; observations of 
a peak in th e brightness of d ebris discs around A stars at ages of 5 
to 10 Myr dCurrie et al . 2008) can be explained in this context but a 
more detailed examination indicates that these models still require 
an inner hole to reproduce the observations. The inner region could 
have been c leared by a planet that is then stirring the outer regions 
of the disc dKennedv & Wvattl20i0h . 

Resolved imaging of young debris discs is therefore important 
in discriminating between discs being perturbed by a fully formed 
planet, such as in the pericentre glow model applied to HR4796 
dWvattetalJI 19991) and disc s in which planet formation may still 
be ongoing such as r\ Tel dSmithetalj [20091) where the disc is 
modelled using the delayed stirring models of iKenvon & Bromlevl 
d2004l) . 
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Many young stars are too distant to study their discs in detail, 
but young moving groups in this age range, such as the BPMG 
(12Myr) provide a relatively nearby, young population in which to 
assess the origin and dive rsity of debris disc emission at this age 
dZuckerman & Sondl2004) . 

HD191089 (HIP 9927 3) is an F5V star at a distance of 
54pc dPerrvman et al.l 1 1997h that was identified as a candidate 
debris disc by [Mannings & Barlowj dl998t) based on IRAS pho- 
tometry. The age of this source is somewhat uncertain and es- 
timates in the li terature vary wi d ely. I sochrone fitting suggests 
old er ages, wit h Nordstrom et alj d2004) giving an age of 3 Gyr 
and IChen et al ] J2006h finding 1.6 Gyr. However, when consider- 
ing evidence of youth such as X-Ray emission and Lithium abun- 
dance data am ongst other techniques, the age indicated is signifi- 
cantly younger; Zuckerman & Songl (l2004) find lOOMyr. However, 
iMoor et al] d2006h suggest that the source is a possible member of 
the /3 Pictoris Moving Group (BPMG) due to HD191089's proper 
motion and high rotation rate, giving it an age of 12 Myr. This mem- 
be rship appears to be confir med from further proper motion studies 
by |Lepine & Simonl (2009). The BPMG is a rich source of debris 
discs with ~23 per cent iRebull et alJl2008h of members harbour- 
ing a disc. Here, we adopt an age of 12Myr when interpreting this 



Using the Spitzer fluxes for this source ^Hillenbran d et al.l 

2008) a nd other sun-like stars , an analysis similar to that pre- 
sented in lSmith & Wvattl ( l2010h for A star debris discs showed that 
HD191089 was the most favourable target for imaging debris discs 
with deep mid-IR observations. 

In §2 we present 18.3(J.m imaging of HD191089 taken with 
the T-ReCS instrument on Gemini South. In §3 we analyse these 
observations and show that we have resolved the debris disc. In 
§4 we confront these observations with models to determine the 
disc parameters and in §5 discuss the implications of the inferred 
structure for the status of planet formation in this system, b oth in 
the context of the Peri centre Glow mod el (Wyatt et "all ll999h and 
the self-stirred model dKenvon & B romley 20041). 



2 OBSERVATIONS 

These observations were taken under proposal GS-2009A-28 using 
T-ReCS on Gemini South with filter Qa (A c = 18.3p.rn, AA = 
1.51pm). The pixel scale of the T-ReCS instrument is 0".09 and the 
total field of view of is 29"x22". These observations were taken in 
parallel chop-nod mode with a chop throw of 15" and chop position 
angle of 55° (East of North). The data were taken over three nights 
(9th June 2009, 12th July 2009, 13th July 2009). The total on source 
time was 7200s. Obse rvations of s t andard star HD 189831 (spectral 
type K5III) listed in ICohen et alj dl999l) . were made to calibrate 
the photometry and also to monitor the PSF of the observations. 
The order of the observations and the on source integration times 
are shown in TableQ] 

These data for HD191089 were taken in sets of four 360s long 
integrations with 6 complete nod cycles ( ABBA). The data we re re- 
duced using custom routines described in ISmith et alj d2008h . The 
data reduction involved determining a gain map for each observa- 
tion using the mean values of each frame to construct a map of pixel 
responsitivity. The on-source pixels were masked during this pro- 
cess, making this equivalent to a sky flat field frame. A DC offset 
was then determined by calculating the mean pixel values in ev- 
ery row and every column, again masking pixels where there was 
source emission present. This was then subtracted from the final 



image to ensure a flat background. Pixels which showed high or 
low gain in comparison with the median response throughout the 
observation were masked off. In order to avoid errors in co-adding 
the data which could arise from misalignment of the images, we fit- 
ted a Gaussian with a sub-pixel centroid to accurately determine the 
centre of the image and so the position of the star. The re-binning 
was done using bilinear interpolation across the array. 

The standard star was observed at a similar airmass to 
HD191089. The calibration levels were compared with the airmass 
for the standard star observations and no correlation was found, so 
no extinction correction was applied to the calibration factors. 

In order to calibrate the images an average of the calibration 
factors determined from standard star observati ons was used. The 
standard star, H P 189831 (F ia = 3.44Jy from ICohen et alJI 19991: 
lAdelmanlk oOl) was observed twice per group of observations (See 
Table Q] for detail of Observations Groups). A coadd of the two 
standards from each group was used to calibrate the observations 
of HD191089 in the same group. The average calibration error for 
the five groups of observations was 6 per cent. The groups had cal- 
ibration factors that varied by 3 per cent, 6 per cent, 12 per cent, 
6 per cent and 4 per cent. These values were used to calculated a 
calibration error for the fluxes. These centred, flux calibrated im- 
ages from each group were then used to produce final images of 
HD 19 1089 and the standard star. 



3 RESULTS 

The final co-added images for HD 19 1089 and the standard star 
HD 18983 1 for PSF reference are shown in FigureQ] The solid line 
contours shown in the images are 25 per cent, 50 per cent and 75 per 
cent of the peak for both HD 19 1089 and the standard. The dashed 
contour represents 10 per cent of the peak emission. The peak value 
for HD191089 in the final coadded QA image is 344mJy/arcsec 2 
and the peak value for the standard is 27648 mjy/arcsec 2 . 

Photometry was performed using a 1".0 radius circular aper- 
ture centred on the star. The stellar flux in the Qa filter is expected to 
be 44mJy from the K band flux (see Section 4.2). Our photometry 
yields a flux of 163±22 mjy including both calibration and photo- 
metric errors for the final co-added image of HD 19 1089 shown in 
FigureQjwith an S/N of 38 (not including calibration errors). Thus 
the photometry gives an excess of 1 19±22 mjy in a 1" radius aper- 
ture centred on the star. The IRS spectrum of lChen et alj2006l after 
subtraction of the stellar flux gives an excess of 115±3mJy. The 
T-ReCs photometry presented here therefore agrees with the IRS 
results. It does not appear that the IRS spectrum includes any emis- 
sion centred on the star outside the 1" aperture used here for pho- 
tometry, as might happen in the larger IRS beam (extracted along a 
3.7" slit). The total error on the photometric measurements in Ta- 
bleQjconsists of both the calibration error and the statistical noise. 
The statistical noise was determined using an annulus with an inner 
radius of 3" and an outer radius of 4" centred on the star, resulting 
in an error averaged over all the groups of l.lmjy/arcsec 2 . 

There is clear ellipticity in the image of HD 19 1089 (centre of 
Figure [Q and a 2D Gaussian fit to HD 19 1089 gives the elliptic- 
ity of the image at the 50 per cent contour as 0.104±0.015 with 
the major axis at position angle of 80° ± 10° East of North. The 
same ellipticity is seen on all three nights. In comparison the ellip- 
ticity of the standard was 0.024 ± 0.007 with the major axis at a 
position angle of 55° ± 2° East of North. As all the observations 
were performed with no on-sky rotation and chop-nod performed 
at 55° East of North, the ellipticity in the PSF could be due to chop 
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Table 1. Observations taken under proposal GS-2009A-Q-28 in order. Note that the on source time is half of the total integration time listed in the table. Fluxes 
are for a 1" radius aperture centered on the star. The group indicates the standard-science-standard observing pattern used. One standard was used in groups 2 
and 3 due to an omitted calibration observation. Each science observation consists of 4 sub integrations. 



Date 


Object 


Group 


Name 


Filter 


Integration Time (s) 


Calibrated Flux (mJy) 


09/06/09 


HD 189831 


1 


Stdl 


QA 


120 


3441 


09/06/09 


HD191089 


1 


Iml 


QA 


1440 


158 ± 21 


09/06/09 


HD 189831 


1 


Std2 


QA 


120 


3441 


12/07/09 


HD 189831 


2 


Std3 


QA 


1440 


3441 


12/07/09 


HD191089 


2 


Im2 


QA 


1440 


167±22 


12/07/09 


HD 189831 


2/3 


Std4 


QA 


120 


3441 


12/07/09 


HD191089 


3 


Im3 


QA 


1440 


166 ± 23 


12/07/09 


HD 189831 


3 


Std5 


QA 


120 


3441 


12/07/09 


HD 189831 


4 


Std6 


QA 


120 


3441 


12/07/09 


HD191089 


4 


Im4 


QA 


1440 


145±24 


12/07/09 


HD 189831 


4 


Std7 


QA 


120 


3441 


13/07/09 


HD 189831 


5 


Std8 


QA 


120 


3441 


13/07/09 


HD191890 


5 


Im5 


QA 


1440 


174±31 


13/07/09 


HD 189831 


5 


Std9 


QA 


120 


3441 



Figure 1. Left to Right: The final co-added images of the standard star (left) and HD 191089 (centre) and a residual emission image shown after subtracting 
an azimuthally symmetric PSF created from the standard star image scaled to the expected photospheric emission from HD 19 1089 (44mJy) from the science 
image (centre). Note that the distinct ellipticity seen in the image of HD191089 is not seen in the standard star. Contours on the left and central plots are at 
25, 50 and 75 per cent of the peak (solid lines) and 10 per cent of the peak (dashed line). The contours on the residual emission image are at 25 per cent, 50 
per cent and 75 per cent of the peak. Orientation of the images is North up, East left. The peak value for HD191089 in QA is 344 mJy/arcsec 2 and the peak 
value for the standard is 27648 mJy/arcsec 2 . The peak value for the residuals is 55mJy/arcsec 2 . The residual image has been smoothed by convolving with a 
Gaussian with a width of 4 pixels. The colour scale is linear from to the maximum pixel value in each image. 
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smearing, but the extension seen in the HD19089 images is at a 
very different position angle and so is unlikely to be an artifact of 
the chopping process. 

The right image in Figure[T]is the residual image of HD191089 
after the subtraction of the standard star image scaled to the flux 
expected from the photosphere of HD 19 1089 (44mJy) followed by 
smoothing by a Gaussian of width 4 pixels. Scaling the subtraction 
to the expected stellar flux means that we are not removing any of 
the disc emission. We also tried subtracting the standard star image 
scaled to the peak of the HD191089 image, to see if there was any 
additional unresolved flux present at the peak, and found that this 
results in a subtraction of 49mJy, indicating that there is very little 
unresolved disc flux contributing to the peak. A ring is seen in the 
residual image, consistent with a disc that is at some intermediate 
inclination between edge on and face on. The position angle of the 
peaks (and hence the major axis of the ring) was found by finding 



the angle between the peaks in the residual image and was found 
to be at 80° ± 10° East of North, consistent with the orientation of 
the ellipticity found for the science observations before subtraction 
of the PSF. 

In order to assess the robustness of this putative disc, and to 
ensure that it is not an artifact of the PSF, we examined the tem- 
poral variability of the FWHM of the science data and assessed 
the significance of the extension. We also looked for any depen- 
dence of the FWHM on time or airmass. We took linecuts at po- 
sition angles corresponding to the major and minor axes of the 
fitted ellipticity in the science images (PA of 80° and 170°) for 
each group of HD191089 observations (see Table [T] for observa- 
tion groups) and for PSF standards from each observational group. 
The linecuts were smoothed by a Gaussian of width 3 pixels and 
were then fitted with a Gaussian profile to determine the FWHM. 
No dependence on airmass or time was found, and the median 
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FWHM values and standard errors were found to be: 8.81±0.2 pix- 
els at 80°for the standard, 9.95±0.3 pixels at 80° for HD191089, 
8.10±0.1 pixels at 170° for the standard, and 9.46±0.4 pixels for 
HD191089 at 170°. There is naturally a large degree of scatter in 
the intensity profiles, due to variations in observing conditions, but 
HD191089 shows extension when compared to the standard at 80° 
and there is also some evidence of extension along the minor axis 
of the image at 170°. This is shown graphically in Figure [2] which 
shows the linecuts through each integration (of 2880s, see table 
[]} of HD191089 (diamonds) and corresponding standard star inte- 
grations (plus signs). HD191089 also has broader wings than the 
standard star image at 80° . 

The residual image of HD 19 1089 shown in Figure [T] also 
shows tentative evidence for a brightness asymmetry between the 
two sides of the disc. The ratio of the flux in a 0.5" radius aper- 
ture, centred on the E peak (found by fitting a Gaussian at sub- 
pixel resolution) and the flux in an identical aperture centred at the 
same radial location on the W side was = 0.80. This aper- 
ture size was chosen to maximise the asymmetry observed. An az- 
imuthally symmetric PSF was used for the subtraction shown in 
Fig.[U so this apparent brightness difference is not simply due to 
non-axisymmetric features in the PSF. This was also checked by 
looking for any systematic azimuthal variations in the PSF, and by 
redoing the subtraction using each of the 10 standards to quantify 
the effect of PSF variability on the measured asymmetry. It was 
found that variability in the PSF FWHM and asymmetries change 
the flux ratio by ±0.02, or < 0.2cr significance. The final stacked 
PSF image used had a FWHM of 8.82 pixels at 80° and 8. 10 pixels 
at 170°. 

The significance of this putative asymmetry was assessed us- 
ing a Monte Carlo method. We created 10 000 axisymmetric disc 
models (see Section 4.1), convolved these models with the axisym- 
metric PSF and added Poissonian noise consistent with the ob- 
served levels in the images (~0. lmjy/pixel). This does not account 
for the effects of any spatially correlated noise, such as coherent 
variation in the PSF over the course of the observations, but when 
we examined the PSF from each group of observations (see Table 
[TJ we found no systematic asymmetries. This modelling provides 
a first guess at the chance of observing such an asymmetry sim- 
ply due to uncertainty from Poissonian noise in the images. The 
average observed asymmetry between the two sides of the image 
in these noisy asymmetric models was 12 per cent. The brightness 
asymmetry observed is then significant at the 1.8<r level. 



4 MODELLING 

4.1 Axisymmetric Disc Model 

To determine constraints on the radial distribution of the emission 
seen in the resolved 18|J.m image we first considered the simplest 
possible model of the disc structure. This model is composed of a 
single axisymmetric disc component, defined by four free parame- 
ters: disc radius (r), width (dr), inclination (i) and surface density 
distribution E, which is assumed to have the form E oc r 1 . The 
disc opening angle (which sets the disc height) is assumed to be 5° 
but this parameter is unconstrained by the modelling process. The 
flux from annuli in the disc at different radii was determined assum- 
ing a grain composition and size distribution that were constrained 
using the emission spectrum as outlined in section 4.2. Since the 
SED fitting process fixes the total flux from the disc model to be 
~115mJy this did not have a significant impact on the derived spa- 
tial distribution. 



Figure 2. The profile of the line cuts through the total images at 80° E 
of N (top) and 170°E of N (bottom). The median FWHM measurements 
taken over all sub-integrations from Gaussian fits to the profiles at the two 
angle are indicated for reference. The averaged profiles from all the line 
cuts for both HD191089 and the standard star are also shown. The pro- 
file for HD191089 (dashed line) at 80° is much broader in the wings than 
that of the standard star (solid line). There is also some evidence of a less 
pronounced extension at 170° as seen in the lower image. Note that mul- 
tiple intensity profiles for the science target (diamonds) and the standard 
(crosses) are shown here. These are taken for each group as indicated in 
table 1. 




Radial offset (arcsec) 



In order to constrain the dust location, as well as the width 
of the disc, the inclination and the surface density profile, a grid 
of models was run. The model images were convolved with an az- 
imuthally symmetric PSF and compared with the observed residu- 
als using the images as well as the linecuts both in the direction of 
extension (PA 80°) and perpendicular to the extension (PA 170°). 
The use of both linecuts allowed us to constrain the inclination si- 
multaneously with the radial morphology. The reduced chi-squared 
of the fit (Xr = ( oi "'~" 1 ° d ) = i where v is the number of free pa- 
rameters and 1 represents a good fit) to each of these three pieces of 
observational data were calculated and were then combined linearly 
with equal weight to come to a final mean reduced chi-squared 

x 2 +x 2 +x 2 
(Xr camiint . d ~ r ''" e8 ° r '" ei7 ° iz±L) w hich was then min- 
imised across the grid of models. The ranges of model parameters 
tested were R in :10 to 50AU (1 AU intervals), R out : 50 to 100AU 
(5 AU intervals), inclination: 10° to 60° (5° intervals, where 0° is 
edge on), surface density index 7: -0.5 to -3.0 (0.5 intervals). The 
values for gamma were chosen to cover possibilities such as the 
surface density distribution expected from grains being blown out 
of the system by radiation pressure (7=- 1.0), and that of the Min- 
imum Mass Solar Nebular (MMSN) (7=- 1.5). The best fit model 
was found to be a ring at 28±1 to 90±5AU, with an inclination 
of 35±5° from edge on, and a surface density profile 7=-1.5±0.5. 
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This best fit model had X? M d = 1-30 with X?iine 80 = l-34, 
Xr!inei 70 = l-40, Xrres = 1.19. The best fit model residual image 
is shown on the left of Figure [3]and can be compared directly with 
the observed disc (Fig[T]right) which has the same colour scale. 

PSF variation can have an important impact on the observa- 
tional residuals, and we used an azimuthally symmetric average 
PSF. This was motivated by the lack of effects of temporal and spa- 
tial variation in the PSF (see Section 3). The modelling process for 
the axisymmetric disc was repeated using the PSF with the largest 
FWHM (Std 8) and the smallest FWHM (Std 6) to find that PSF 
variation had no effect on the modelling results within the quoted 
errors. 



4.2 SED Models 

The emission spectrum of HD191089 is shown in Figure [4] The 
stellar photospheric emission has been modelle d using a Ku rucz 
model atmosphere with a temperature of 6440K jKuruczl 19791) . fit- 
ted to the 2MASS K band photometry implying a stellar luminosity 
of 3.2Z/0 at the inferred distance of 53.5 pc. Photospheric flux con- 
tributions have been subtracted from the fluxes shown at > 5 urn 
which thus represent the excess emission from the disc. This spec- 
trum is approximately that of a blackbody fit with a temperature of 



99 ± 9K and a fr actional luminosity (/ = ±^i) of 1.5 x 10~ 3 
dChen et al]|2006h . 

To calculate the emission from grains in the model they were 
assumed to have a size distribution with n(D) oc D -3 ' 5 , where D 
is the grain diameter, whic h is the standard solution for a theoreti- 
cal collisional cascade ("see lDohn anvi 1969), that is truncated at a 
minimum and maximum grain size. For the assumed composition a 
grid of models was calculated with a range of minimum sizes 0.5- 
IODm where Dm is the largest grain size removed from the system 
by radiation pressure for which (3=0.5 where j3 = p™ d n , the ra- 
tio between radiation and gravitational forces (Burns et al .11 19791) 
which depends on stellar properties, grain composition and poros- 
ity. The maximum size was fixed at 1cm since larger grains have a 
negligible contribution to flux from a model with this size distribu- 
tion. The grains were assumed to have a silicate core (amorphous 
olivine) and an accreted mantle of organic refractories produced 
by UV photoprocess ing of ice. (as used in lLi & Greenberjl 19971 ; 
Auger eau et al.1 19991) . A range of compositions was also tried, with 
amorphous silicate fractions varying from per cent to 90 per cent 
by volume and with porosities (i.e. vacuum fraction of grain by 
volume) from per cent to 95 per cent. Di electric constants were 
calculated from tabula ted laboratory values dLi & Greenberd 19971 ; 
Augereau et al. 1999) using Maxwell-Garnett effective medium 
theory. The optical properties of the grains were calculated using 
Mie theory, Rayleigh- Gans theory and Geometric Optics in the ap- 
propriate size regimes (Bohren & Huffman 1983). 

This resulted in a grid of 360 models over which \ 2 was 
minimised by fitting to the spectrum and MIPS fluxes (24|im and 
10\im). The fit to the IRS spectrum was calculated using Xirs ~ 



ad) 



where F a t s is the observed flux, F mo d is the model 
flux and N is the number of points in the observed IRS spectrum. 
This normalises Xirs so that when calculating the final \ 2 the fit 
to the spectrum was given equal weight to the fits to the 24|im and 
70vim MIPS fluxes. The best fit model with a reduced \ 2 of 1.14 
has a minimum grain size that is coincident with the blow-out grain 
size ( 1 .46 vim for the composition used) and a composition of amor- 
phous silicate fraction of 10 per cent by volume, and a porosity of 
60 per cent with the rest of the grain composed of organic refrac- 



tory materials, with no ices present. The fractional luminosity is 
1.4 x 10~ 3 . As there are no obvious spectral features in the SED, 
the main constraint on the composition and size distribution comes 
from allowing the appropriate range of temperatures to be present 
given the constraint that the dust is in the region 28-90 AU; black- 
body grains would have to be at a radius of 15 AU to achieve the 
observed temperature. Although our model provides a consistent 
fit to both image and SED, it is not expected that the composition 
has been uniquely constrained by this process. This dust model was 
then used as an input to the model grid and there was no change to 
the best fitting model of the 18|^m images. 

The total mass in the collisional cascade, with the assumed 
size distribution of n(D) oc D -3 ' 5 scales as M to t oc y/D max - 
Thus for Dmax = 1km the total mass would be 13M^ two or- 
ders of magnitude hig her than the current mass of the Kuiper Belt ( 
iBernstein et alj |2004^. but note that we have no information on the 
maximum size of objects in this system. 



4.3 Emptiness of the Inner Hole 

It is possible to fit all excesses in the spectrum with a model that 
only extends from 28AU-90AU. There is no requirement either 
from the image or from the spectrum for emission at <28AU. The 
observed 18(J.m image of HD191089 shows a resolved belt of dust 
from 28-90 AU, with a cleared inner region extending to a radius 
of ~28AU. We compared the photometry of the observations in the 
inner 28AU (0.5" radius aperture) and then repeated the measure- 
ment with the same aperture on the models, where an unresolved 
flux component could be added at the position of the star. Excluding 
calibration uncertainties, the measured flux density within a 0.5" 
radius aperture on the 18u.m image was 48±4mJy of flux, with 
a photospheric contribution of 33mJy (75 per cent of total photo- 
spheric flux, measured in the same aperture using a scaled PSF). 
The same measurement on the model with no stellar or unresolved 
flux component added gives 1 6m Jy. Therefore the maximum unre- 
solved component that can be added to the model that does not 
exceed a 3a upper limit of the observations is 12mly. However 
these measurements do not account for the 12 per cent calibration 
uncertainty in the photometry or the error in photospheric flux due 
to uncertainties in the spectral fitting to the star, which is ±1.2mJy 
(on total photospheric flux). Since a change in calibration would 
also result in a corresponding change in the mode fluxl (to fit the 
peaks at the correct level), it is instructive to consider what would 
have been derived with a revised calibration factor and stellar flux. 
In the extreme situation that the calibration was 36 per cent higher 
(i.e. a 3<r upper limit on the calibration factor) and the stellar flux 
was lower by 3.2mJy (3a lower limit), the flux measured in the 
previously described apertures would have been 65±5mJy with a 
corresponding model flux of 22mJy. Thus the excess would have 
beenl3±5mly and so gives a maximum unresolved component of 
<28mJy for a 3a deviation in both calibration and statistical uncer- 
tainty. 



4.4 Pericentre Glow Model 

The brightness asymmetry se en in Figure [T] is re miniscent of that 
seen in the disc of HR4796 dTelesco et al.l l2000). The 6 per cent 
asymmetry of HR4796 was interpreted as a consequence of peri- 
centre glow, where the secular gravitational perturbations of a 
planet on an eccentric orbit imposes a forced eccentricity on the 
planetesimal belt, causing the forced pericentre side of the disc 
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Figure 3. The best fitting axisymmetric disc model for the 18|^m emission (see section 4.1). The model residual image (leftO can be compared directly with 
the observed residuals (Figure[T] right). The difference between the observed residuals and the model residuals (right) indicates an over subtraction near the 
disc ansae and a slight under subtraction at the peak. See section |5~Tl for details. 
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Figure 4. Spectral Energy Distribution (SED) of HD191089. The photosphere of HD191089 is fitted with a Kurucz model profile (M sta r = 1.4M , 
Lstar = 3.21/0, T a tar = 6440-R") scaled to the 2 MASS K ban d flux and shown with a solid line. Observations at > 5u.m are shown after the subtraction of 
the photospheric contributions: the IRS spectrum of Chen et al. 2006 (solid line from 5 — 30)J.m), MIPS fluxes (Carpenter et aL 2009) (triangles) and IRAS 
colour-corrected excesses (asterixes) and Gemini TReCs fluxes in 0.5" radius aperture (square), and sub-mm/mm upper limits (crosses) from SHARC and 
CARMA. The excess measurements are first fitted with a 99K modified blackbody (dashed line), for comparison to a realistic grain fit (solid line) described in 
section 4.2 
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to be closer to the star and hence hotter and brighter (Wy att et"al] 
jl999i) . Here we applied the same model to the HD 19 1089 disc. 
The disc is defined by the inner and outer semi-major axis of the 
disc, &min and St max with a proper eccentricity e p =0 (so a. m i„ 
and a ma x are equivalent to i m in and t max in the axisymmetric 
model), the power law index for the distribution of semi-major axes 
IPG = 7+1 where 7 is the surface density index from the axisym- 
metric model, and a forced eccentricity e / . The best fit model has 
parameters a m j T1 — 28, a ma x— 90, 

7pg = 2.5 and e/=0.12 and is 

shown in Figure[5] 

This addition of an extra parameter (the forced eccentricity) 



improves the Xr from 1.31 to 1.15, but in general adding 

an extra parameter will improve the fit of any model. However, 
given the tentative (1.8cr) nature of the asymmetry, is adding this 
extra parameter to the model justified? We assess the justification 
for this extra parameter through use of the Bayesian Information 
Criterion (BIC): 

BIC = Nln(Xr combintd ) + klnN 

where N is the number of data points, k is the number of 
free parameters and Xr aombined ls me minimum combined reduced 
\ 2 for the model. The BIC considers the fit of the model but 
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Figure 5. The best fitting disc from the pericentre glow model which is scaled and shown in the same way as the HD 19 1089 residuals (Figure [T} for 
comparison. The disc parameters are the same as those for the best-fitting axisymmetric disc shown in Figure [5] (28-90AU, i=35°, Eocr 7 ) Right: linecut 
through the observed and model images at a PA of 80° (Observed: solid line, Model: dotted line) 
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penalises the model for extra parameters(see IWahhai et al]|2005l : 
lLiddlell2007t) . Two models were considered: the axisymmetric disc 
model(see section |4~H and the pericentre glow model (see section 
14.4b . The BIC value for the axisymmetric disc is 19.3 and for the 
pericentre glow model is is 19.2. There is very little difference be- 
tween the BIC values, indicating that there is no preference be- 
tween the two models. 



5 DISCUSSION 

5.1 Implications for the Planetary System 

There are 3 different possible interpretations of the state of planet 
formation in the HD191089 system based on the resolved image: i) 
The inner cavity and brightness asymmetry are due to the formation 
of planets interior to the inner edge of the resolved disc (J5.1.1t , ii) 
The disc is in the process of forming Pluto-siz ed bodies as in the 
delayed stirring models of I Ken von & Bromley] d2005h J5.1.2t or Hi) 
The disc is a remnant of the protoplanetary disc d5.1.3t . 

5.1.1 Planet in the Central Cavity 

The 18|J.m image of HD191089 shows a resolved belt of dust from 
28-90 AU, inclined at 35° from edge on with a cleared inner region 
extending to a radius of ~28AU. The brightness asymmetry seen in 
the HD 19 1089 disc could be evidence for a planet on an eccentric 
orbit (see Section |4~4l (. The pericentre glow model only places con- 
straints on the forced eccentricity imposed on the debris disc by the 
planet's secular perturbations and this cannot be used to constrain 
its mass and semi-major axis. The only constraint is that the lower 
the mass of the planet the longer the secular precession timescales. 

Although it is not known it seems likely that should such a 
planet exist, the same planet that is causing the as ymmetry is caus- 
ing the inner edge, as proposed bv lOuillerj 1200(3) . The inner edge 
of the resolved dust belt at 28AU could therefore mark the outer 
edge of the planetary system, especially if a giant planet is sculpting 
the inner edge of the disc , as in the case of Fomalhaut ( Kal as et al.l 
l2008tlChiang et all2009» which has a dust belt with a well defined 
inner edge and a brightness asymmetry. 

The steepness of the inner edge of the disc also suggests trun- 
cation by a planet. For a disc that extends inward to ~5AU, like 



the discs seen in iKenvon & Broml ey (2004), a surface density of 
Ear 6 from 5-28AU is required to satisfy the observed 3a upper 
limit of 27mJy (see section 4.3) on emission from this region. This 
is steeper t han surface density p rofile of Eoc r 7 ^ 3 predicted in the 
models of IKenvon & Bromlevl ((2004). For a surface density pro- 
file like those of the delayed stirring models to be consistent with 
the observations, then the disc would have to be truncated at 24 AU. 
Assuming the disc is being truncated by a planet at just inside 28AU 
and that the disc is currently being stirred, then this constraint on 
the secular timescale means a minimum mass f or the perturber can 
be derived using Eq. 42 from l Wvatt etail d 1 999» which gives a min- 
imum planet mass of 2.6Me for a system age of 12 Myr. 



5.1.2 Delayed Stirring Model 

Another possible interpretation of these observations is that the re- 
solved belt corresponds to a region of a more extended disc that has 
recently formed Plu t o size d objects. In these self-stirred models of 
IKenvon & Bromlevl d2005f) the evolution of an extended planetes- 
imal belt is followed, allowing both the velocity and size distri- 
butions to evolve self-consistently. A destructive collisional cas- 
cade is only ignited in the disc when planets forming in the disc 
reach Pluto-size. The timescale for the formation of these Pluto 
sized objects depends strongly on their distance from the star and 
on the surface density of the disc so the outer disc regions are 
stirred at later times. We performed furthe r modelling b ased on 
the prescriptions for the evolution given in Wvatt (2008) and re- 
fined in Kennedy & Wyatt (2010) which provides an empirical fit 
to the self-stirred models of Kenyon & Bromley (2005). This model 
consists of a planetesimal belt with radii r min and i ma x with sur- 
face density scaled to the MMSN divided into annuli at different 
radii which have a suppressed dust production rate until the time at 
which Pluto sized bodies have formed, stirring the annulus into a 
destructive collisional cascade, increasing the production of small 
dust and causing it to brighten. The dust emission is calculated as- 
suming emission from dust grains with a size distribution and com- 
position consistent with the SED and image modelling in previous 
sections and a ratio of dust area to planetesimal mass appropriate 
for a collisional cascade size distribution. To account for deviations 
from this theoretical size distribution a scaling factor r\ was applied 
to the surface brightness of the whole disc. However, under these 
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constraints from the SED and with the stirring time set to the age 
of HD 19 1089 (12Myr) the model provided a poor fit to the obser- 
vations. This is primarily due to the shape of the surface density 
distribution produced.To produce enough emission in the outer re- 
gions of the disc (i.e out to ~90AU) there has to be a region of 
gradually rising surface density (with a power law form Eoc r 7//3 
) in the inner region of the disc (see Kennedy & Wyatt 2010 Fig. 
2). For reasons similar to those discussed in 5. 1 . 1 this overproduces 
emission at <28 AU (0.5"), where the 3a upper limit on the flux in 
this region is 27mJy. Therefore the delayed stirring model cannot 
fit the sharp inner edge to the disc. However this can be achieved 
if there is an additional mechanism truncating the disc. For exam- 
ple a disc with an inner radius of 24 AU, surface density profile 
E oc r 7/3 between 24 and 34 AU, and E oc r" 1 ' 5 from 34 to 90 
AU can fit the observations. This gives a flux of 26mJy in the inner 
28AU (0.5") of the disc, consistent with the 3<r upper limit. We con- 
clude that some mechanism must be truncating the disc, and that it 
is not possible to determine from these observations alone whether 
planet formation is ongoing in some or all of the outer disc. 



5.1.3 Protoplanetary Disc Remnant 

A further possibility is that this disc is a remnant of the proto- 
planetary disc and the dust is confined in a ring due to interac- 
tions with gas that remains in the system from the protop lanetary 
phase dTakeuchi & Artvmowiczll200ll ; iBesla & Wull2007t) . These 
models are too gas rich for HD 19 1089 as a massive gas disc 
around a star of this age would be unusual as most gas dis cs are 
thought to dissipate in <10Myr (e.g Holle nbach et alj|2005l ). Gas 
has been detected around 4 A stars known to have debris discs (/3 



Pic: iBrandeker et alj|2004l 49 C eti: bent et aljl2005l , HD141569: 
Jonkheid et alj|2006l HD32297: lRedfieldll2007f) and its effect on 



Krivov et al.l 1 120091) concl uded that gas does no t af- 



dust is debated: 

feet the dust distribution in /3 Pic, though IBesla & WuJ J2007I) dis- 
cussed how gas may be undetected and still have an effect on the 
dust dynamics. Deep gas observations towards HD 19 1089 would 
be needed to rule out this effect. 



5.2 Placing HD191089 into context: Comparison with F stars 
of the BPMG 

The BPMG (Beta Pictoris Moving Group) is a collection of young 
(-~12Myr) stars with common proper motion, indicating that they 
formed in the same cluster which is now dispersing. Their veloc- 
ities can be traced back to this common point of origin. The age 
of the Moving Group can be assessed using the M star members, 
and the group is then assumed to be coeval. Table [2] gives a cen- 
sus of the currently known F star members of the Pic M oving 
Group. Membership was taken from lLepine & Simonl (2009). who 
assessed the membership of the moving group by examining stellar 
proper motions. By studying the fraction of these stars that have in- 
frared excesses we can place HD191089 into context with its coeval 
moving group. Of the nine known F star BPMG members, 6 have 
infrared excesses, 3 of which ha ve resolved debris disc s : HD15115 
jKalas et alj|2007l) , HD181327 ^Schneider et"ai1l2006l : Ichen et al.l 
l2008h andHD191089 (this work). 

For the 6 stars with discs SED fitting was used to estimate 
radius and fractional luminosity, which are given in Table 2, and 
range from 10 to 150 AU, and from 3 x 10" 5 to 3 x 10" 3 . 
HD 19 1089 is the second brightest and has the second smallest inner 
radius. It is the most similar to HD181327 which is slightly brighter 



and larger. The similarity is also apparent in ima ges. HD181327 
has b een imaged in both scattered light using HST {Schneider et alj 
2006) and in mid-IR thermal emission dChen et all2008t This disc 
shows an extended icy "Kuiper belt" at 86AU with a width of 36AU 
dSchneider et alj2006l) . compared to a belt at 55AU with a width of 
50AU for HD191089. Both these discs have cleared inner regions 
(HD181327 ~68 AU, HD191089 ~28AU) and both show evidence 
for a brightness asymmetry between the two sides of the disc in the 
mid-IR (HD181327: flux ratio of 1.4±0.1 between the Northern 
and Southern arms HD191089: flux ratio of 0.8±0.12 between the 
E and W sides of the disc). 

The other resolved disc HD15115 is a very large, highly asym- 
metric disc, dubbed the "Blue Needle" due to the elon gation of the 
West side of the disc to >550AU (East side ~ 310AU) jKalas et all 
120071) . Such an extreme asymmetry could be due to a stellar fly-by, 
and HIP 12545 (also a member of the BP MG) is a po ssible candi- 
date due to its near on-sky position jKalas et al] |2007). This asym- 
metry could also be caused by interaction with a nearby ISM cloud, 
that could b e 'sandblasting' and so stripping mater ial from one side 
of the disc jDebes et al.ll2009l: iManess et al]|2009h . There are not 
many obvious similarities between HD 19 1089 and HD15115: the 
unusual morphology of HD15115 appears to originate through in- 
teraction with an external perturber rather than due to any plan- 
ets in the system and it appears much more extended in images 
than HD191089. Vega's disc displays a similarly extended shape 
to HD151 15, with a high populatio n of sma ll grains in the process 
of being blown out of the system dSu et al. 2005). However, this 
system can still be explained in the context of a steady state col- 
lisional cascade when the difference in stellar luminosity between 
the pole and equat or due to Vega's fast rotation is taken into account 
dMUller et ai]|2010l) . 

The remaining three stars are all binaries, and one of these, 
HD199143, may show some evidence for an infrared exces s 
due to a circumsecondary disc ( Jayawardhana & Brandeker 200 lh . 
but this may simply be due to uncer tainties in the photospheric 
modelling of the binary companion 1 van denAnckeretaljbood: 
IChauvin et al]|2002l : iKaisler et alj|2004h . The general lack of ex- 
cess seen for the F stars of the BPMG that are also bina- 
ries could be due to perturbations from the secondary truncating 
any possible discs. HD199143 has a seperation of ^ 53-105AU 
jjavawardhana & Brand eker 200j]; IChauvin et alj|2003) and there 
is observational evidence that binary companions in the range 50- 
100AU are expected to disrupt discs (Jensen et al. 1994. The other 
two binarie s have quite w ide predicted separatio ns ( HD 14082: 
~530 AU ISong et all [20031 . HD29391: ~800AU iFeigelson et al.l 
2006), and unless on very eccentric orbits the secondaries are un- 
likely to preclude the presence of discs around the primaries in 
these systems. It is possible that the binary has evolved or that there 
is something about the formation mechanism of binary systems that 
suppresses formation or detectability of a debris disc, but the small 
number statistics for this sample means that this is still inconclu- 



6 CONCLUSIONS 

We have presented Qa band Gemini T-ReCs imaging of the 12 Myr 
old F star HD191089. The emission at 18.3(im is shown to be sig- 
nificantly extended compared with a point source. This image rep- 
resents the first resolution of dust emission around this star. These 
observations confirm the interpretation of the SED as a debris disc 
with a single component resolved at 28-90 AU. There is little emis- 
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Table 2. Census of F star members of the BPMG. The dust radius is calculated by fitting a black body to the 24)^m and 70|xm fluxes, or from previous literature 
black body fits to the SED. 



Star 


Spectral 


Mass (M ) 


Distance 


Luminosity 


L Dust 
L» 


Blackbody Dust 


Refs 


Notes 




Type 




(pc) 


(£©) 




Radius 
[Resolved] (AU) 






HD203 


F3V 


1.3 


39.1 


2.7 


1 X 10~ 4 


10 


1,2,3,4 


1 


HD 14082A 


F5V 


1.23 


39.42 


1.88 


N/A 


N/A 


1,2,3,4 


2 (Binary) 


HD 15115 


F2V 


1.39 


44.8 


2.94 


4.9 x 10" 4 


34 [-300-500] 


5,6,7 


3 (Resolved disc) 


HD 29391 


FOV 


1.27 


29.9 


5.85 


N/A 


N/A 


1 


4 (Binary) 


HD 35850 


F7V 


1.19 


26.8 


1.78 


3 x 10~ 5 


150 


1,8,9,10,11 


5 


HD 164249 


F5V 


1.32 


46.8 


2.64 


5.9 x 10~ 4 


20 


1 


6 


HD 181327 


F5/6V 


1.3 


51 


3.1 


3xl0~ 3 


20 [86] 


7,1,6,12,13 


7 (Resolved disc) 


HD 191089 


F5V 


1.4 


54 


3.2 


1.4x10" 3 


15 [55] 


1,7,8,14 




HD 199143 


F8V 


1.2 


48 


2.45 


7 


7 


15,16,17,18 


8 



References: 



Rebull et al. 2008 


Car] 


renter et alj2008l 




Lepine & Simonl 



Carpenter et al. 2009 
Kalas et al.|2007l 
Moor et al. 2006 



7 Chen et al. 200i 



8. n^occatagliata et al l 
120091 

9. ? 



10. 
11. 

12. 

13. 
14. 



Spangl er et alj200ll 



Silverstone et al 
Chen et alj2008 




Mosh ir et al][l993 
IChen et all200d 
15 Jayawardhana & Brandekei 2001 
16 lKaisler et alj2004l 
17. I van den Ancker et alj2000l 
18. lchauvin et al.l2002l 



Notes: 

1. Strong Spitzer 24[im excess (-F24 



60F>94) ICarpenter et alJl200S 



2. Binary companion HD 14082B (G5V, 530Au separation). ISong et*aT 
20031 



3. Resolved Asymmetric Debris disc " Blue Needle". Possibly perturbed 
by HIP12545 flvbv. lKalas et alj2007l 

4. Binary companion (-800AU). Blended with S pitzer. 

5. Dust based on Spitzer 70|xm excess. Isilverstone etalj gOO^ 

6. Dust based on Spitzer 24|xm and 70um exc ess. iReb ull et al. l2008l 

7. Resolved dis c at -86AU (-20AU width) I Schneider et alj 120061 : 
IChen et alj2008l 

8^ Bin ary compani o n (M 2, 86AU [javawardhana & B randekerl 

l200ll; iKaisler et al.1 120041). Possible c ir cumse co ndary disc 
ivan den Ancker et al l l200(t IChauvin et all I2002L IKaisler et all 
l2004h . 
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sion from the inner 28AU of the disc, indicating the presence of an 
inner cavity. 

The disc also shows a tentative brightness asymmetry of 
■p^- = 0.80 ± 0.12 with a significance of 1.8a. This asymme- 
try is consistent with a scenario in which an perturbing body on an 
eccentric orbit imposes a forced eccentricity of 0.12 on the plan- 
etesimal belt through secular perturbations, causing the centre of 
the ring to become offset from the star which replicates the bright- 
ness asymmetry observed. This interpretation can be tested by fur- 
ther observations seeking to confirm the brightness asymmetry and 
measure the predicted offset of the star from the centre of the ring. 

To predict the surface brigh tness of the HD191089 disc in 
scattered light we used S = £^ dWeinberger et al.1 19991) where S 
is the surface brightness (mjy/arcsec 2 ), F is the stellar flux (mjy), 
r is the optical depth, uj is the albedo and (j> is the angular dis- 
tance of the disc from the star. Using our best fit model parameters 
at NICMOS wavelengths (A c =l.l, AA=0.59) we predict a surface 
brightness of 0.024-0. 19mJy/arcsec 2 assuming isotropic scattering 
and an albedo of 0.1-0. 8 the range predicted for Kuiper belt objects 
JStansberrv et al.ll2008h . Although this is detectable with ACS and 
NICMOS resolving the disc would be difficult as the peak surface 
brightness is expected at a radius of 1", which will be close to the 
inner working edge of the coronograph for NICMOS observations, 
and inside it for ACS. 

Three models for the central hole were considered. The in- 
terpretation that fits most neatly with the model for the bright- 
ness asymmetry is that the inner region has been cleared of plan- 
etisimals by a planetary system. In this case high resolution im- 
ages of the system could resemble Fomalhaut, including a planet 
close to the inner edge of the belt. We found that the sharpness of 
the inner edge was incomp atible with the delayed stirring models 
dKenvon & Bromlevl [2004) in which the disc is in the process of 
forming Pluto sized objects causing the disc to brighten at the radii 
where they have just formed, although it is possible that planet for- 
mation of the type envisioned by these models is ongoing as long 
as some other process is truncating the disc at 25AU. A third possi- 
bility is that the planetesimal belt is somehow confined due to dy- 
namical interactions with gas in the system. However, there are no 
observational constraints on the presence of gas in the HD 19 1089 
system. When compared to other members of the BPMG, the most 
obvious similarities in terms of fractional luminosity and radius in- 
ferred from blackbody fitting are with HD181327, whi ch has also 
been resolved at mid-IR wavelengths dChen et alj 2008) and shows 
evidence for an inner clearing and asymmetry similar to those seen 
in HD191089. 
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